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Introduction 

Neural  coincidence  detection  is  essential  in  sound  localization,  which  requires  the 
computation  of  interaural  time  differences  (ITDs)  (for  frequencies  below  a  few  kHz).  This 
is  performed  by  binaural  cells  in  the  avian  Nucleus  Laminaris  (NL),  and  its  mammalian 
homologue,  the  Medial  Superior  Olive  (MSO). 

A  “coincidence  detector”  neuron  should  fire  when  inputs  from  two  independent  neural 
sources  coincide  (or  almost  coincide),  but  not  when  two  inputs  from  the  same  neural  source 
(almost)  coincide.  A  neuron  that  sums  its  inputs  linearly  would  not  be  able  to  distinguish 
between  these  two  scenarios.  Segregating  the  inputs  on  separate  dendrites  should  avoid  this 
problem:  post- synaptic  depolarization  arising  from  a  synaptic  event  is  reduced  if  the 
dendrite  has  already  been  previously  partially  depolarized.  This  idea  has  been  used  by 
Agmon-Snir,  Carr  and  Rinzel  to  model  bipolar  dendrites  as  interaural  coincidence  detectors 
in  the  avian  auditory  brainstem.  This  project  is  a  more  biophysical  model,  using  the 
program  NEURON,  to  explore  the  robustness  of  this  mechanism  and  look  for  new 
mechanisms. 


In  addition,  the  active  dendritic  potassium  channels  (both  low  and  high  voltage  activated), 
are  incorporated.  The  high  voltage  activated  channels  (Kvl.2)  vastly  increase  the  cell’s 
ability  to  phase  lock  at  frequencies  >  1  kHz. 
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Model  Description 

The  model  emulates  a  single  neuron  with  an  axon,  soma,  and  a  variable  number  of  dendrites, 
each  section  with  a  variable  number  of  equipotential  compartments.  All  geometric, 
electrical,  and  channel  parameters  are  adjustable,  as  are  the  number  of  synapses/dendrite  (~ 
30),  the  synaptic  locations,  and  the  distribution  of  synaptic  locations.  Channel  types  include 
potassium  (high  and  low  voltage  activated  [Kvl.l,  1.2],  and  delayed  rectifier),  sodium,  and 
passive.  Active  channels  have  adjustable  time  and  voltage  parameters. 

The  stimulus  is  a  pure  tone  of  adjustable  frequency,  with  variable  binaural  phase  difference 
(or  contralateral  monaural  stimulus  with  variable  ipsilateral  spontaneous  activity).  More 
complex  stimuli  can  be  easily  introduced. 

The  synapses  fire  with  conductance  proportional  to  an  alpha-function,  with  adjustable  time 
constant,  peak  conductance,  and  reversal  potential.  The  synapses  fire  as  individual  Poisson 
processes,  with  probability  rate  given  by  the  half-wave  rectified  pure-tone,  with  adjustable 
amplitude  and  base  spontaneous  firing  rate.  The  fast  Kv3.1  channels  of  the  pre-synaptic 
neurons  are  reflected  in  a  short  synaptic  time  constant. 


The  implementation  uses  the  program  NEURON  and  has  an  straightforward  graphical  user 
interface  for  controlling  the  parameters  and  running  the  model.  NEURON  provides  a  real¬ 
time  display  of  data  and  analysis  including  the  potential  at  various  locations,  the  two  stimuli, 
the  synaptic  firings,  spike  counters,  period  histograms  of  synaptic  firings  and  the  action 
potentials,  and  their  vector  strengths. 
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NEURON  Panels 
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Spatial  Geometry 

This  typical  model  cell  has  2  dendrites, 
each  ~60  pm  long  and  4  pm  in  diameter, 

a  spherical  soma  of  diameter  15  pm,  and  ^ 

an  axon  100  pm  long  and  2  pm  in 

diameter. 


“Space  plots”  show  the  intracellular  potential  at  every  point  in  the  cell 
(here  at  one  instant  in  time). 


Potential  along  Axon,  Soma,  &  Dendrite[0] 


The  potential  along  the  curve: 
down  the  axon,  through  the  soma, 
and  down  along  one  dendrite. 


C  enter  for  Aud Kory 
and  Acoustic  Research 


Potential  up  Denrite[0]  &  down  Dendrite[1] 


The  potential  along  the  curve: 
up  one  dendrite,  through  the  soma, 
and  down  along  the  other  dendrite. 
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Time  Plots 


Typical  time  plots  (continuously  updated).  This  pair  of  cells  each  receives  the  same  stimulus 
probability  distributions  (at  500  Hz).  The  top  cell  receives  its  inputs  binaurally  in-phase,  and 
the  bottom  out-of-phase. 


The  red  curve  tracks  the  intracellular  potential  at  soma/axon  boundary,  the  yellow  at  the 
axon  tip,  and  the  brown  curve  in  one  dendrite.  Under  are  a  pair  of  plots  of  the  presynaptic 
probability  distribution.  The  bottom  8  curves  of  each  graph  show  actual  synaptic  currents 
(note  the  spread  from  the  Poisson  distribution). 
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Period  Stimulus  Histograms 


Spikes 


Ipsilateral  Inputs 


Contralateral  Inputs 


360 


360 


The  same  pair  of  cells  (and  stimulus),  tracked  for  250  ms.  The  vector  strength  (VS)  is  a 
measure  of  phase  locking,  ranging  between  0  and  100%.  The  cell  receiving  in-phase 
stimulus  increased  its  vector  strength,  the  out-of-phase  cell  substantially  reduced  it.  Note 
also  the  number  of  spikes  in  each  case. 
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Parameter  Space  &  Dendritic  Length 


The  dimensionality  of  the  parameter  space  of  this  model  is  high.  Only  a  small 
volume  of  the  space  of  it  is  biologically  relevant,  but  it  is  not  immediately  obvious 
where  the  relevant  subspace  lies. 

The  parameters  fall  into  a  few  general  categories:  those  whose  values  are 
experimentally  known,  those  whose  values  are  not.  Some  values  fall  in  a  particular 
range,  and  different  ranging  parameters  may  be  correlated  or  not.  Some  parameters 
may  be  particularly  relevant  for  certain  species  but  not  for  others  (e.g.  the  barn  owl 
can  detect  ITDs  up  to  8  kHz,  whereas  the  chicken  can  detect  only  up  to  2  kHz). 

To  limit  the  search,  we  usually  compare  the  performance  of  a  pair  of  identical  cells  to 
identical  stimuli,  with  one  set  of  stimuli  in-phase  and  the  other  out-of-phase. 

As  a  first  step,  we  pick  “reasonable”  values  for  all  parameters,  vary  the  dendritic 
length,  compare  firing  rates  and  observe  the  vector  strength  of  the  in-phase  case. 


The  “nominal”  firing 
rate  of  -350  spikes/s  is 
a  typical  best  firing  rate 
for  an  NL  cell  in  chick. 


In-Phase  Rate 
Out-Phase  Rate 

In-Phase/ 
Out-phase 
Ratio  [%] 

In-Phase 
Vector 
Strength  [%] 
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Stimulus  Frequency 


Next  hold  fixed  the  dendritic  length  at  which  the  in-phase  firing  rate  matches  the 
“nominal  rate,”  and  vary  the  frequency. 


This  cell 
responds  best 
to  frequencies 
between  300 
and  1200  Hz. 


In-Phase  Rate 
Out-Phase  Rate 

In-Phase/ 
Out-phase 
Ratio  [%] 

In-Phase 
Vector 
Strength  [%] 


In  the  avian  NL,  the  dendritic  length  varies  (roughly)  inversely  with  the  best  frequency  of 
the  cell.  By  directly  covarying  the  stimulus  frequency  with  the  dendritic  length,  we  can 
attempt  to  capture  a  central  subspace  of  the  entire  parameter  space. 


stimulus  frequency  [Hz] 


With  these 
parameters  this 
class  of  cells 
does  poorly 
below  800  Hz. 
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Canonical  Parameters 


Covarying  the  stimulus  frequency  with  the  dendritic  length  gives  us  a  curve  that  is 
relatively  flat,  but  by  covarying  one  additional  parameter,  the  maximum  dendritic 
conductance  of  the  high  voltage  activated  potassium  channel  (KHvaX  we  can  get  a  rate 
curve  that  matches  the  nominal  rate  quite  closely  (at  least  at  250  Hz  and  above). 
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Synaptic  Distribution 


Typically  the  synapses 
would  be  uniformly 
distributed  along  the 
dendrite.  We  can 
compare  the  base 
response  to  the  response 
when  all  synapses  are 
moved  to  the  center  of 
the  dendrite. 

The  top  figure  is  the  base 
response,  the  bottom 
with  all  synapses  at  the 
center. 

The  effect  is  relatively 
small.  The  In-Phase  rate 
is  moderately  reduced, 
worsening  at  large 
dendrites/low 
frequencies.  The  ratio  of 
In-Phase  rate  to  Out- 
Phase  rate  remains  large 
in  the  same  range. 
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One  might  have  expected  the  cell  to  perform  better  with 
all  the  synapses  together,  to  reduce  jitter  from  the 
differing  travel  times  of  the  signals  to  the  soma,  but  the 
dendrites  are  (electrotonically)  relatively  compact. 
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Synaptic  Location 


We  can  also  change 
the  location  of  the 
concentrated  synapses. 

Here  are  the  results 
from  3  cases:  all 
synapses  located  at  the 
dendrite  base; 
halfway  down  the 
dendrite;  at  the 
dendrite  tip. 

At  the  base,  the 
coincidence  detector 
performs  poorly  for 
large  dendrites/low 
frequencies,  giving 
many  false  positives. 
The  advantage  of 
inputs  adding 
sublinearly  in  the 
dendrites  (which  are 
much  smaller 
compartments  than  the 
soma)  is  lost. 
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Synaptic  Reversal  Potential 


The  benefit  of  sub- 
linear  addition  can  also 
be  seen  by  changing 
the  synaptic  reversal 
potential. 

Here  are  the  results 
from  3  cases:  setting 
the  reversal  potential 
10  mV  too  low, 
correctly,  and  10  mV 
too  high. 

When  the  reversal 
potential  is  artificially 
raised,  the  sub-linear 
addition  of  the  inputs  is 
suppressed,  and  the 
coincidence  detectors 
become  effective  over 
a  smaller  range  of 
frequencies/lengths . 
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Dendriteless  Soma 

This  model  supports  the  claim  that  dendrites  are  aiding  in  neural  computation.  How 
well  does  it  perform  compared  to  the  same  model  without  dendrites? 

In  the  previous  “Synaptic  Location”  panel,  if  the  synapses  are  located  at  the 
dendrite/soma  boundary,  the  performance  is  compromised.  But  the  dendrites  are  still 
present  and  can  still  act  as  current  sinks,  which  can  increase  coincidence  detection 
ability. 

In  this  trial  the  dendrites  have  been  removed. 


4000  2000  1000  500  250  125 
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Firing  rates  rise  dramatically,  but  they  do  so  for  both  the  In-Phase  and  Out-Phase  cases. 
The  ratio  of  the  rates  is  near  1,  indicating  poor  discrimination. 
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Spontaneous  Input 


In  this  model,  the  input 
vector  strength  (VS)  is 
determined  by  the 
spontaneous  activity. 

s  =  0  gives  no  spontaneous 
activity.  The  half-wave 
rectified  inputs  have 
VS  =  1/V2  =  77%. 
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high  frequencies. 
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Number  of  Synapses/Dendrite 


Coincidence 
detection  is  robust 
against  changes  in 
the  number  of 
synapses/dendrite. 
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Number  of  Dendrites 


Real  NL  Neurons  have 
more  dendrites  at 
shorter  lengths/higher 
frequencies. 

In  this  model,  more 
dendrites  may  enhance 
coincidence  detection 
at  high  frequencies, 
but  the  evidence  is 
murky. 
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Real  NL  Neurons 
vary  in  dendritic 
diameter  from  2  to 
4  |Lim. 

The  major  effects 
of  this  variation  are 
changes  in  the 
electrotonic  length 
and  in  the 
transmission 
coefficient  to  the 
soma. 

This  model 
indicates  no  benefit 
to  having  the 
thinner  dendrite. 
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Synaptic  Refractoriness 


stimulus  frequency  [Hz] 


The  refractory  time  of 
the  pre- synaptic 
inputs  may  be 
important  for  the  low 
frequency/long 
dendrite  case. 
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The  NL  contains 
arrays  of  coincidence 
detectors  tuned  to  a 
wide  range  of  phase 
differences  (more  than 
just  0°  and  180°). 

Here  are  responses  to  a 
360°  range  of  phase 
difference. 

The  firing  rate  and 
vector  strength  are 
both  strong  functions 
of  phase  difference, 
but  the  firing  rate  is 
more  sharply  tuned 
(and  the  vector 
strength  is  not  reliable 
at  low  rates). 
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Phase  Locking 


Rate  [spikes/s]  Vector  Strength  [%] 


1587  Hz 


250  Hz 
95  |lin 


3175  Hz 


500  Hz 
57  pm 


0  90  180  270 


198  Hz 
113  pm 


2520  Hz 


397  Hz 
67  pm 


0  90  180  270 

157  Hz 
135  pm 


2000  Hz 


315  Hz 
80  pm 


0  90  180  270 

125  Hz 
160  pm 


Interaural  Phase  [°] 


Institute  for  Systems  Research 
University  of  Maryland 


Monaural  Input 

Real  NL  neurons  also  respond  to  periodic  stimuli  presented  only  in  one  ear. 
Presumably  this  is  due  to  coincidence  detection  between  the  monaural  stimulus  one 
side  and  the  spontaneous  activity  on  the  other  side. 

The  model  cell  locks  well  to  this  monaural  input,  with  a  vector  strength  of  near  1  for  a 
broad  range  of  inputs,  though  its  firing  rate  is  well  below  the  nominal  best  rate  for 
binaural  input. 
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Conclusions 


•  The  model  has  parameter  ranges  that  give  behavior  corresponding  to  the 

behavior  of  real  NL  neurons. 

•  The  dendrites  aid  in  the  ability  of  the  cell  to  perform  coincidence  detection, 

especially  because  of  sublinear  addition  and  dendritic  current  sinks. 

•  The  high  voltage  activation  potassium  channels  are  important  for  coincidence 

detection  at  high  frequencies  (though  roughly  independent  of  the  kinetics). 

•  Coincidence  detection  is  robust  against  the  number  of  incoming  synapses. 

•  Regions  in  the  parameter  space  have  been  identified  for  further  investigation. 

•  The  model  predicts  that  vector  strength  is  very  robust  (at  fast  firing  rates),  but 

is  not  as  sharply  tuned  as  firing  rate  as  a  function  of  phase  difference. 

•  The  model  locks  well  to  monaural  stimulus. 

•  NEURON  is  an  extremely  effective  environment  for  implementing  a  detailed 

biophysical  model. 
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Abstract 

The  coincidence  detector  neurons  in  the  auditory  brainstem  detect 
interaural  time  differences  (ITDs).  Inspired  by  a  minimal 
biophysical  model  [Agmon-Snir  et  al.  Nature  393,  268-272  (1998)], 
we  constructed  a  model  using  the  NEURON  environment.  ITD 
coding  was  improved  when  the  inputs  from  both  ears  were  located 
on  the  bipolar  dendrites  and  segregated  over  having  both  inputs  on 
the  soma.  Thus  the  model  behaves  both  like  the  in  vivo 
coincidence  detectors  and  the  minimal  model.  The  model  has 
enabled  us  to  explore  features  of  the  coincidence  detector  neurons 
unexplained  by  the  simple  biophysical  model,  including  the  effect 
of  synapse  location  and  multiple  dendrites. 
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